Mass freezing in growing neutrino quintessence 
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Growing neutrino quintessence solves the coincidence problem for dark energy by a growing cos- 
mological value of the neutrino mass which emerges from a cosmon-neutrino interaction stronger 
than gravity. The cosmon-mediated attraction between neutrinos induces the formation of large 
scale neutrino lumps in a recent cosmological epoch. We argue that the non-linearities in the cos- 
mon field equations stop the further increase of the neutrino mass within sufficiently dense and large 
lumps. As a result, we find the neutrino induced gravitational potential to be substantially reduced 
when compared to linear extrapolations. We furthermore demonstrate that inside a lump the possi- 
ble time variation of fundamental constants is much smaller than their cosmological evolution. This 
feature may reconcile current geophysical bounds with claimed cosmological variations of the fine 
structure constant. 

PACS numbers: 98.80.-k,98.80.Jk 
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I. INTRODUCTION 

The theoretical and observational implications of a 
coupling between a light quintessence field or cosnion and 
other matte r sp ecies have been widely explored in the 
literature [ll-[l5|. However, both cosmological tests and 
laboratory measurements challenge many of these scenar- 
ios and lead to strong bounds on the allowed strength of 
the coupling. This motivates the search for mechanisms 
that tend to naturally suppress the observable traces of 
these interactions which could be manifested, for exam- 
ple, in deviations from the gravitational force law or in 
the time-dependence of fundamental constants. 

One possible suppression mechanism relies on the den- 
sity dependence of fundamental couplings due to the den- 
sity dependence of the value of a light scalar field [l6| . 
If the mass of the scalar field is larger than the inverse 
size of a dense object, then the local value of the field 
typically adapts to the local density of the object rather 
than to the cosmological evolution. For the example of 
the chameleon mechanism [T7l - [20| the mass of the scalar 
field itself increases with density such that for sufficiently 
large coupling and dense objects the scalar field value de- 
couples from the cosmological evolution. The scalar field, 
called a chameleon field under these circumstances, is al- 
lowed to interact with all forms of matter. The strength 
of the coupling can be substantially larger than that of 
gravity. Strong enough non-linearities are a general in- 
gredient for this type of mechanism to work. 

Growing neutrino quintessence [111, [2^ is also charac- 
terized by important non-linearities. The coupling be- 
tween the cosmon and neutrinos is substantially larger 
than the gravitational coupling and results in the forma- 
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tion of very large neutrino structures at late times j23| . 
In this work, we focus on this scenario and discuss the 
implications for the gravitational potential of such neu- 
trino lumps and on the variation of the electromagnetic 
fine structure constant in the presence of a field coupling 
to electromagnetism. 

We find a strong backreaction effect which substan- 
tially reduces the neutrino mass within such neutrino 
lumps as compared to its cosmological value. In conse- 
quence, this effect leads to a substantial reduction of the 
neutrino-lump induced gravitational potential. Inside a 
neutrino lump the possible time variation of fundamen- 
tal constants is strongly suppressed as compared to their 
cosmological time evolution. 

Strong non-linear features are present in this scenario 
despite the fact that the density dependent mass of the 
scalar field typically remains smaller than the inverse size 
of the object and that the chameleon mechanism does not 
operate. The main ingredient of the non-linear behavior 
is the non-linear dependence of the local neutrino mass 
on the local value of the cosmon field. This combines 
with a non-linear cosmon potential. 

We assume that a very rapid increase of a local fluc- 
tuation in the neutrino number density ends after some 
virialization process in a quasi-static lump containing a 
given number of neutrinos. For subsequent time we take 
a static neutrino number density n^(f) which is a func- 
tion of the distance to the center of the lump. We then 
solve the field equation for the cosmon field (f> taking into 
account the cosmological time variation of (f) outside the 
lump. The field equation for the local cosmon field (j) 
has a source term proportional to I3m^(<j))n^{r), where 
j3 is the cosmon-neutrino coupling and may depend on 
(j). The value of the neutrino mass m^{(f)) depends on (jj 
and we assume, for simplicity, an equal mass for all three 
neutrino species. In the scalar field equation, the source 
term may be partly balanced by gradient terms and the 
derivative of the cosmon potential, dV{(j))/d(j). This re- 
sults in an almost static value of the cosmon field which 
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depends, however, on the neutrino density n^{r). 

A key ingredient of the growing neutrino quintessence 
scenario is the coupUng of the cosmon to the trace of 
the stress energy tensor of neutrinos, {p^, — 3p,y), where 
Pi, and denote the neutrino energy density and pres- 
sure, respectively. For most of the cosmological evolu- 
tion neutrinos are relativistic and almost massless, thus 

~ Pi,/3. Consequently, in this regime, the coupling is 
ineffective and both components (cosmon and neutrinos) 
evolve independently. The relativistic neutrinos perme- 
ate the universe homogeneously, only subject to cosmic 
expansion as in "standard cosmology" . Furthermore, the 
cosmon tracks the background evolution following an at- 
tractor solution on trajectories characteristic of the pres- 
ence of an exponential field potential 0, [23 - l27| . How- 
ever, as soon as the neutrinos turn non-relativistic, the 
interaction switches on as ~ 0. Since the abundance of 
cosmic neutrinos is predicted by standard cosmology, the 
moment in time when neutrinos become non-relativistic 
is determined by the present value of the mass as well 
as the strength of the cosmon-neutrino coupling and the 
cosmon potential. Typically, in the growing neutrino sce- 
nario, the cosmon-neutrino coupling exceeds the strength 
of the gravitational interaction by two to three orders of 
magnitude. The neutrino mass grows more rapidly for 
larger values of the coupling, but starts to increase later 
in time such that neutrinos become non-relativistic at 
redshift znr ^5 — 10. 

As the neutrino mass grows, in turn, the neutrinos gain 
influence on the dynamics of the scalar field. Eventually 
the cosmon is forced out of the attractor and the early 
cosmological solution crosses over to an almost constant 
asymptotic field value at late times. Consequently, the 
energy density corresponding to the value of the expo- 
nential cosmon potential at the transition, acts similarly 
to a cosmological constant and drives cosmic accelera- 
tion. In the growing neutrino scenario, the correspond- 
ing dark energy scale is directly related to the particle 
physics scale of the present neutrino mass [111, . This 
is a possible solution of the "why now?" problem of dark 
energy. 

Only after neutrinos have turned non-relativistic, the 
cosmon becomes the mediator of a "fifth force" between 
neutrinos. The strength of this interaction is propor- 
tional to the square of the coupling. The range of this 
force is of the order of the Hubble length, since the time 
evolving mass of the cosmon roughly equals the Hub- 
ble parameter. This is in contrast to adiabatic models 
of mass varying neutrinos (MaVaNs), which assume a 
scalar field with a mass much larger than the Hubble 
scale [H-Ull . The mass of the cosmon emerges naturally 
from possible explanations of an exponential potential 
stemming from an asymptotically vanishing dilatation 
anomaly [33j . 

Due to the strong long-ranged force, at z < znr, fiuc- 
tuations in the non-relativistic neutrinos grow rapidly on 
all scales larger than the free-streaming length. Recent 
studies have widely explored the evolution of the grow- 



ing neutrino perturbations both in the linear [34| - l36l | and 
in the non- linear regime ^1], finding a substantial clus- 
tering of neutrinos up to supercluster scales and beyond. 
Typically, neutrino perturbations with sizes larger than 
>10 Mpc become nonlinear at a redshift of z w 1 — 2 
(SJ-jSg. Around this redshift the universe is thus pop- 
ulated by virialized neutrino lumps bound by the "fifth 
force" and by gravity. Present findings indicate that the 
neutrino profile within these lumps is consistent with a 
highly-concentrated "Navarro-Frenk- White" (NFW) dis- 
tribution. 

Noting that the impact of these overdensities on the 
local scalar field dynamics is amplified by the large cou- 
pling, one can ask whether neutrino lumps provide a new 
source for highly non-linear effects. In this paper we ex- 
plore this question and analyze the resulting theoretical 
and observational consequences. To this end, we first de- 
termine the spatial cosmon profile within individual viri- 
alized neutrino lumps embedded within a time-evolving 
homogeneous neutrino background. 

By considering example cases, we find the non-linear 
structure of the field equation to have the stunning effect 
of decoupling the value of the cosmon within a lump from 
its large-scale cosmological behavior. This implies a sub- 
stantial reduction of the neutrino mass within neutrino 
lumps. To probe the dependence of these effects on the 
properties of the neutrino lumps, we investigate a fam- 
ily of neutrino distribution functions ni,{r). For a given 
n^{r) we numerically solve the scalar field equation, with 
the boundary condition that outside the object the cos- 
mon asymptotes to the cosmological value. This allows 
us to determine the average neutrino mass in the interior 
of various families of neutrino lumps. The average neu- 
trino mass turns out to decrease with increasing number 
of neutrinos in the lump. This implies a substantial re- 
duction of the gravitational potential $,y of large scale 
neutrino lumps compared to linear extrapolations which 
are based on the cosmological neutrino mass. These find- 
ings reflect the importance of backreaction effects in the 
growing neutrino scenario. Furthermore, they bear sig- 
nificant consequences for the effect of neutrino lumps on 
the large scale structure and the cosmic neutrino back- 
ground in the form of the Integrated Sachs- Wolfe effect 

Finally, one may wonder whether the described non- 
linear effects in growing neutrino cosmologies might im- 
ply new interpretations of cosmological observations. For 
example, we consider the implications of our findings 
for the variation of the fine structure constant a which 
would result from a coupling of the cosmon to electro- 
magnetism. We demonstrate that within large neutrino 
lumps, the non-linear effects suppress the time variation 
of the fine structure constant, while allowing for a strong 
variation in the background. This may reconcile geophys- 
ical bounds with claimed cosmological observations of a 
time variation of a. 

This paper is organized as follows. After setting the 
stage for our analysis, in Sec. |lT] we set up the general 
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equations governing the cosmon dynamics inside and out- 
side of a neutrino lump. In Sec. IIIII we present our nu- 
merical results for the spatial cosmon profile and the neu- 
trino mass and comment on the difference in field behav- 
ior compared to linear approximations. In Sec. IIVI we 
present our results for the average neutrino mass for var- 
ious families of lumps and comment on the implications 
for their gravitational potential. We also discuss possi- 
ble backreaction effects on the background cosmological 
value of the scalar field and the growth of large scale 
neutrino perturbations. In Sec. |V]we explore the conse- 
quences of our finding for the variation of the fine struc- 
ture constant in the presence of a scalar field coupled to 
electromagnetism. Finally, in Sec I VII we summarize our 
results and conclude. In a brief appendix, we derive re- 
lations between the parameters characterizing the NFW 
neutrino distribution function. 



II. SETTING THE STAGE 

In this work, it is our aim to estimate the spatial cos- 
mon profile at late times for various families of neutrino 
lumps, formed at a redshift 2:nl — 1 — 2. To this aim we 
consider the perturbation to the field generated by a sin- 
gle virialized neutrino lump embedded in a time-evolving 
neutrino background of uniform density. This approxi- 
mation is justified in the limit that the distance between 
lumps is much larger than their virial radius Tvir. 

As we will show in Sec. Illli for a given static neutrino 
distribution in the interior and a given time-evolving 
background density well outside, the properties of a lump 
and its impact on the field can essentially be quantified by 
a single parameter, namely, the total number of neutrinos 
Nu it contains. In this section, we will present the main 
ingredients of the growing neutrino scenario pTI . [23 [ and 
describe the background evolution of the cosmon. Fur- 
thermore, we will set up the field equation for the spatial 
dependence of the cosmon within the approach outlined 
above. 



A. The growing neutrino scenario 

We take both the scalar field potential V{4)) as well as 
the neutrino mass mi,{(j}) to depend exponentially on the 
scalar field value. 



V(0) = exp 



M 



m,y{(f>) = rh exp ( —jS 



M 



(1) 



(2) 



where M = {SttGn)'^^^ denotes the reduced Planck 
mass, TO is a constant that fixes the neutrino mass scale 
and a > and /3 < are free model parameters. To com- 
ply with observational bounds on the amount of early 
dark energy, the parameter specifying the exponential 



potential is constrained to be a > 10 [4^. Note that ac- 
cording to our choice in Eq. the dimensionless cou- 
pling /3 generally defined as 



-M- 



(3) 



is a constant. However, similar effects to the ones we 
encounter for a constant coupling are expected in the case 
of a field dependent coupling, 0((t>) , as has been suggested 
within a particle physics context [22|. 

As in previous works (2ll . [23| . we assume a = 10 
and /3 = —52 for our numerical studies. We choose a 
present cosmological neutrino mass TOj^q = 2.3eV consis- 
tent with current cosmological and experimental bounds. 
The value of rh in Eq. ^ is then dictated by the precise 
timing of the transition from matter to cosmon domina- 
tion according to the observed dark energy abundance 
17do « 0.73. 



B. Time evolution of the cosmon in the 
cosmological background 

The dynamics of the cosmological background field 
is described by the Klein Gordon equation. 



3p.(0)), (4) 



where dots denote the derivative with respect to confor- 
mal time t and Ti — a/a is the corresponding Hubble 
parameter. Furthermore, pi, and pi, denote the neutrino 
energy density and pressure, respectively. One observes 
that the neutrino source term on the right hand side of 
Eq. (lU is strongly suppressed as long as neutrinos are rel- 
ativistic, since ~ Pi//3. However, when they turn non- 
relativistic, Pi, ~ 0, the source term Pp„/M becomes of 
dynamical importance. For constant /?, the energy den- 
sities of the cosmon and of the neutrinos evolve asymp- 
totically proportional to each other. Averaging over os- 
cillations in time, the cosmological value (jj evolves as 



(z) = 00 



3M 

a — (3 



M , .... , 
— In (1 



ln(l + z) , 

a//3) 



(5) 



with z being the cosmic redshift and pdo = V{4') + Pi' 
V{(t)) denoting the dark energy density. 



C. Quasi-static cosmon field in a neutrino lump 

Initial neutrino perturbations can evolve into virialized 
lumps with a highly concentrated NFW type neutrino 
distribution. Due to the neutrino-cosmon coupling, the 
spatial profile of the scalar field (j>{r) resulting from a 
given static spherically symmetric distribution n^(r) of 
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non-relativistic neutrinos is determined from the follow- 
ing field equation: 



r d(h M 



where 



(6) 



(7) 



Here and in what follows a prime denotes a derivative 
with respect to the radial coordinate r. Furthermore, 
the last equality applies for non-relativistic neutrinos, 
with pi, and denoting the neutrino energy density and 
number density, respectively. 

In growing neutrino quintessence, the source term in 
the spatial field equation depends of the field value in 
contrast to other classes of coupled dark energy scenar- 
ios [41]. The reason is that m((/)) in Eq. ^ cannot be 
linearized as (3(j)/M ;» 1 owing to the large coupling and 
the large field value. This issue will be discussed in detail 
in Sec. HIl 

The solution of the field equation ^ has to obey the 
following boundary conditions 



4>'{r)\ 



and 4'{r)\r- 



(8) 



The first condition follows from requiring the field value 
to be finite and well defined at r = and the second 
implies that the field must converge to its cosmological 
background value </> in the limit of r going to infinity. The 
possible dependence of the solution on redshift is induced 
by the redshift dependence of the background field value 

We assume here, implicitly, that the neutrino number 
distribution n^{r) for a virialized lump has decoupled 
from the cosmological evolution and therefore is time- 
independent. In addition to the static lump distribu- 
tion n„{r), the cosmon also couples to the cosmological 
background neutrino density n^{z) which is much smaller 
than ny{r) within the lump. This coupling to the back- 
ground, therefore, only plays a role at the border or out- 
side the lump. 



D. Neutrino distribution 

The spherically averaged radial distribution of neutri- 
nos within virialized lumps is typically described_by a 
highly-concentrated Navarro-Frenk- White profile [231] 



?^NFw(?') 



(r/rs) (1 -I- r/rs 



(9) 



where is the normalization of the number density 
which can be different for different lumps. The char- 
acteristic scale separates the inner part r <^ with 
riNFW from the outer region, where riNFW oc . 



The number of neutrinos contained within a radius R can 
be written as 

N^{R)^4tt drr'^n^Fwir) = N^FiR), (10) 
Jo 



with 



F{R) = \n(l + - 



R 



R + Ts 



(11) 

(12) 



Because the increase of Ni^{R) for large R ^ Vg is only 
logarithmic. 



N^{R)c^N,ln{ — 



(13) 



we can use as a reasonable parametrization for the 
number of neutrinos in the lump. More exactly, ap- 
proximately gives the number of neutrinos within a ra- 
dius lOrs. 

This implies that instead of Us and one may use 
and Vs for the parametrization of the profile n^, (r) = 
n-MFw{f)- We typically compare with the total num- 
ber of neutrinos within our horizon. 
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where n,yo denotes the present cosmological number den- 
sity of neutrinos and Hq the present Hubble parameter. 
Thus, for a NFW profile, we will characterize the neu- 
trino lumps by the two numbers Af*/iVtot and r^. 

It is also possible that the logarithmic increase of 
Ni,{R) does not extend to arbitrary large R. This hap- 
pens if a neutrino lump competes for the neutrinos in 
neighboring lumps. One may imagine that the infall of 
further neutrinos is stopped at some characteristic time 
and subsequently all available neutrinos concentrate in a 
finite volume of size characterized by a virial radius rvir. 
We account for this possibility by considering a second 
type of profile 



nstcp(r) 



r > r„ 



(15) 



In what follows we will refer to these two possibilities 
Eqs. (HI) and ^ as the "pure NFW" and the "step 
NFW", respectively. We will see that the pure NFW 
(which corresponds to rvir — >■ oo) and the step NFW 
yield qualitatively similar results. 

The step NFW profile ([T5|) implies a well specified 
number of neutrinos corresponding to the lump 



N^ = ATr dr n^Fw{r) F{c) , (16) 

Jo 



with 



F(c) = ln(l + c) 



l + c 



lnc-1, (17) 



5 



where c = Vvir/rs and the last equahty holds for c ^ 1. 
A neutrino lump is in this case characterized by three 
parameters: A^*/7Vtot, fs and the concentration c. For a 
finite Tvir or c we may also replace iV^/TVtot by Ni^/Ntot- 
Note that for each class of lumps with given N^, and 
rvir the maximal value for the number of lumps A// within 
the horizon is 



(18) 



This corresponds to the limit when all neutrinos in the 
visible universe are captured in identical lumps. 

In order to satisfy the second of the boundary condi- 
tions in Eq. (jS]) we must add to the lump profile a cos- 
mological background number density ni,{z) such that 



with 



n^{z, r) = riNFwl'") ^'(^vir ~ r) + ni,{z) , (19) 



n,{z)^n,o{l + zf , (20) 



where n^o corresponds to the present value and we re- 
call that Tvir — > oo for the pure NFW profile. Of course 
Eq. (I^Hl) only remains a reasonable estimate for the neu- 
trino number density provided most neutrinos are not yet 
concentrated in lumps. Once most of the neutrinos are 
in lumps, the time evolution of the cosmological average 
value of the quintessence field (/>(z) is more complicated 
and backreaction effects ought to be taken into account. 
For our purposes Eq. ([20]) is sufficient, since we are only 
interested in the effects of a given time evolution 4'{z) 
on the cosmon profile inside a lump and not in the com- 
putation of the evolution of 0(z) itself. In general, the 
redshift dependence of 4){z) can be made consistent by an 
approximate formal choice of n^{z) in Eq. (|19p . The de- 
tailed redshift dependence of cf) and n^{z) is unimportant 
for the findings of this work. 



III. MASS FREEZING 

In this section we demonstrate that the neutrino mass 
inside of sufficiently dense and large neutrino lumps does 
not follow the cosmological evolution. For this purpose, 
we numerically solve the cosmon field equation, Eq. ([S]), 
for spherically symmetric static neutrino distributions Q 
and ([T5|) . The boundary condition 0'(r = 0) = ensures 
the regularity of the solution. The second boundary con- 
dition is given by the asymptotic value of (f> for large 
r — > oo, (/)(r — oo) = (see Eq. (|5])). The value 
must equal the cosmological average value of the cosmon 
field (or more precisely, the average over regions with- 
out lumps). Thus we compute the solutions for different 
"boundary values" 0. The time or redshift dependence of 
the cosmon field in the lump is then induced by inserting 
the cosmological value of (/>(-z). This procedure assumes 
that the adaptation of the cosmon field inside the lump 
to the cosmological value occurs on a sufficiently small 
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FIG. 1: Scalar field profile 4>{r,z = 0) in the radial direction 
for a pure NFW neutrino profile (Eq. (|9|) and for a step NFW 
profile (Eq. QSll). For both cases = 4xlO"^iVtot. We show 
curves for = l.SMpc and = 2.6Mpc. For the step NFW 
profile we set c = 30. 



time scale compared to the inverse Hubble parameter, 
and that the use of a static solution provides a good ap- 
proximation. In practice we impose the boundary value 
(j){z) for some very large value of r ^ Ts where the deriva- 
tive of the potential and the source term in the field equa- 
tion nearly cancel each other. 

Let us start by discussing the approximation of linear 
field equations, where the field dependence of the neu- 
trino mass m,y{(l>) in Eq. (6) is neglected and dV/d(j> « 0. 
In this case, the general solution can be written as 
(j}{r) = (j) + S4>{r), where 5<j){r) vanishes for r 00 and 
is independent of 0. A cosmological change of (j) would 
then result in a displacement of the cosmon field for the 
whole lump configuration. Such a linear approximation 
would be valid for uncoupled quintessence or sufficiently 
small /? - in this case the cosmon field acts as a "cosmo- 
logical clock" . If the variation of fundamental couplings 
(such as the fine structure constant) can be written as a 
linear dependence on the displacement of 0, then, their 
evolution inside the lump follows the cosmological evolu- 
tion [H-Iil. 

In our case of strong coupling, (3^1, the linear 
approximation is not justified as the ^-dependence of 
rrii, oc exp(— /30/M) plays an important role. We demon- 
strate this here for the case of a constant /3. For a 0- 
dependent /?(</>) the effect may even be more pronounced. 
In Fig. [T]we illustrate the radial dependence of 4>{r) for 
the pure and step NFW neutrino profiles for a lump with 
A^* = 4 X lO^^A^tot, evaluated for ^0 = (i>{z — 0) corre- 
sponding to the present cosmological value of the field. 
The values of the parameters will be motivated at the 
end of this section. We find a qualitatively similar be- 
havior 4>{r) for the two profiles. In the outer region of 
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Pure NFW 
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FIG. 2: Neutrino mass profile m^{r)/fh,^o in the radial di- 
rection for redshifts 2 = and z = 0.14 for the pure NFW 
profile. 



Step NFW 




FIG. 3: Neutrino mass profile m^(r)/m^o in the radial di- 
rection for redshifts z = and z = 0.14 for the step NFW 
profile. 



the lump, the value of is not important. As one might 
have expected, the deviations of 0(r) from the cosmolog- 
ical value (j) occur for larger r for the pure NFW profile. 
In the inner region of the lump mainly matters (for a 
given Nf), while the value of r^ir only plays a minor role. 
Though the typical deviations of (/)(r) from the cosmo- 
logical value are only half a percent, their impact on the 
value of the neutrino mass is of considerable importance 
due to the large value of the coupling (3. 

For the various profiles of n„{r), we found the value 
of nil, in the interior to be roughly three to four orders 
of magnitude smaller than the present cosmological one, 
rhi^Q. This significant reduction owes to the exponential 
dependence mi,{r)/mi, = exp(— /3(0(r) — ^)/M), which 
makes the neutrino mass very sensitive to the decrease 
of the field value in the interior of a lump. The impli- 
cations of this result for the gravitational potential of a 
neutrino lump as well as for the cosmological evolution 
of the cosmon will be discussed in Sec. IIVI 

These findings suggest that the neutrino mass inside 
the lump does not follow the cosmological evolution. We 
demonstrate this in Figs. [2] and [31 where we evaluate the 
neutrino mass for the same neutrino profiles of Fig. [T] 
and for two different redshifts z = and z = 0.14. We 
compare the neutrino mass as a function of the distance 
from the center of the lump with the present cosmological 
value of the neutrino mass, rh^Q. For large r the neutrino 
mass follows the cosmological evolution independently of 
the properties of the lump. In contrast, for small r the 
change of the neutrino mass is very substantially reduced. 
The curves for the same profile ni,{r) but with different 
4>iz) are nearly on top of each other. Inside the lump, the 
value of m,j(r) depends essentially only on the parame- 
ters characterizing the neutrino number density profile of 
the lump, mainly on N,, and to a smaller extent, on r^. 



This demonstrates very directly the "mass freezing" ef- 
fect. The time evolution of the neutrino mass inside the 
lump is frozen at its value at some redshift characteristic 
for the formation of the lump. It depends on the neutrino 
density of the lump, but not on the cosmological evolu- 
tion. This decoupling of the local properties from the 
cosmological evolution is very similar to the decoupling 
of the gravitational potential of a galaxy from the cosmo- 
logical expansion. We will discuss in Sec. |V] that this de- 
coupling can have important consequences for a possible 
observation of time variation of fundamental constants. 

We demonstrate the decoupling explicitly in Fig. SI 
where we show the difference of the local field between 
z = and z = 0.14, A0(r) = = 0, r) - = 0.14, r) 
for the neutrino number densities described above. The 
figure illustrates that in the considered example cases the 
radial field profile in the interior changes up to three or- 
ders of magnitude less than for large r. It can be verified 
that a larger iV* enhances this effect. 

We conclude that non-linear effects turn out to be very 
strong and it is interesting to compare the behavior of 
the solution within the lump with a "local density ap- 
proximation" where the gradient terms in the scalar field 
equation are neglected. In this case, the cosmon field 
locally follows the neutrino number density according to. 



(21) 



or equivalently 



/3 



M 



^(^-^"'/^'m.on.lr). (22) 



Here we choose a constant 00 a-nd assume that l^(0o) 
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FIG. 4: Time change of the cosmon profile A(^(r) = (f>{z = 
0, r) — </)(2 = 0.14, r) in the radial direction for four neutrino 
number density profiles. 



equals the present dark energy density, 

V{4>o) - 0.7pc ~ 10-120^4, 



(23) 



where pc denotes the present critical energy density. We 
compare in Figs.[5]and|6]the numerical solution of Eq. ^ 
with the local density approximation (j22[) for various neu- 
trino number profiles (r) . We observe that though the 
local density approximation and the numerical solution 
are convergent for large r, they are fairly different for 
small r. This is a reflection of the fact that the cosmon 
mass is not large enough for the chameleon mechanism to 
operate in our growing neutrino quintessence model [4l| . 

In other words, the boundary conditions ([8]) prevent 
the field to sit at the minimum of the effective potential, 
i^min, but instead at larger values. The field cannot be 
at smaller values 0(r = 0) < (j^min, because this would 
lead to a runaway solution (if wc think of r as the time 
variable) and the boundary condition at infinity would 
not be reached. For the local approximation to work, the 
mass of the cosmon would have to be much larger than 
the inverse size of the lump. 

We finally comment on our choice of for a given 
value of A^*. Even though the dependence of our numeri- 
cal results on Ts is moderate, as can be seen by comparing 
Figs. 2 and 3, it is substantial enough to warrant some 
estimate of this parameter. For a given N^,^ the radius Ts 
determines the density in the core of the lump. Realizing 
that nNFw('"s) = "-5/4, one has nj^Fw{rs) = N^/lGirr^. 
For r ^ one finds fT.NFw('') — nsTs/r ~ N^./Airr'^r. 
Thus smaller enhances the core density. This qualita- 
tive feature also applies to lumps with a finite size rvir, 
if we keep N^, fixed (both for fixed c or fixed rvir). 

Some insight on orders of magnitude can be gained 
from the early stages of the formation of individual lumps 
which has been analyzed in Ref. [1^. Given that the 
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FIG. 5: Comparison of the local density approximation (thin 
lines) with the numerical solution (thick lines) for the scalar 
field profile for a pure NFW neutrino profile. 
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lines) with the numerical solution (thick lines) for the scalar 
field profile for a step NFW neutrino profile. 



cosmon-neutrino interaction is ineffective during the rel- 
ativistic neutrino regime, the initial power spectrum for 
the linear neutrino perturbations is the same as in stan- 
dard cosmology. Once the neutrinos turn non-relativistic, 
the cosmon-neutrino interaction becomes effective and, 
as a result of the large coupling, the neutrino pertur- 
bations grow rapidly on all scales larger than the free- 
streaming length. 

After a first phase, where an initial neutrino fluctua- 
tion accompanies the universe's expansion, its potential 
energy becomes important with respect to the kinetic 
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FIG. 7: The number of neutrinos N^, within a sphere of ra- 
dius R at two epochs during the coUapse of a lump. The 
hnes represent the corresponding fits to a pure NFW-profile. 
The data on were obtained using results of a dedicated 
numerical analysis. 



fitted this particular case, which corresponds to a final 
gravitational radius i?f = 72.1Mpc, to the NFW profile 
(cf. Eq. Uni)). We define R{ as in Ref. [H, i.e. the dis- 
tance to the center of the lump where the gravitational 
potential is a factor 2 smaller than the potential at the 
center. 

For Zi = 2.044 one finds iV*/iVtot = 1-2 x I0~^ = 
64.2Mpc, while Zf = 2.017 yields N^/Ntot = 5.5 x 10"^ 
Vg = I5.7Mpc. It is clear that at these redshifts the 
neutrinos continue to fall into the structure, resulting 
in an increase of iV* and with z. Nevertheless, this 
gives already an idea of the order of magnitude of these 
quantities for larger lumps. Once the further infall of 
neutrinos is stopped or substantially reduced, e.g. by 
competing lumps, we expect the concentration process of 
the neutrinos to go on for a while until a static virialized 
situation is reached. This concentration process leads to 
a reduction of for almost constant . This motivates 
the values of used in our figures. For lumps with a 
finite radius Tvir we may also relate to the neutrino 
number density contrast A„ = (n)/n^o- This is discussed 
in the appendix. 



energy of the neutrinos it contains. Eventually, the per- 
turbation reaches its maximal size and thereafter, with 
increasing inward velocities, contracts according to the 
laws of free fall until virialization occurs and the collapse 
terminates. 

The Navier-Stokes equations describe the formation 
process of an individual neutrino lump in an expanding 
universe in physical space as a function of time. A nu- 
merical code has been developed 23|, which allows one 
to track the corresponding time evolution of the relative 
neutrino density perturbation Si, = Spy/pv as well as 
the cosmon perturbation Scf) = 4)(r) — cj) from the linear 
regime up to the highly non-linear regime, until shortly 
before the neutrino lump virializes. In this code, there is 
no assumption on the final number density profile of the 
lump. The neutrino number density distribution is re- 
lated to the spherically averaged density contrast S^, and 
field perturbation Scj) at a given redshift Zf by, 

n.(r, Zf) - n,o(l + z,f{l + 5,(r, Zf)) e^^^(-'-f)A'^ . (24) 

Using results of this numerical analysis, we show in Fig. [7] 
the number of neutrinos (R) within a sphere with ra- 
dius R, 



(25) 



Ni,{R) = 47r / drr^n^ir), 



at two moments in time: when the code terminates 
(zf = 2.017) and shortly before (zi = 2.044). (The termi- 
nation of the code due to numerical instabilities can be 
associated with the onset of virialization. At this moment 
the potential and the kinetic energy of the neutrinos in 
the lump have reached roughly equal values.) We have 



IV. GRAVITATIONAL POTENTIAL OF 
NEUTRINO LUMPS 



The strength of the non-linear effects depends on the 
number of neutrinos in a lump as we will demonstrate in 
this section. For this purpose we consider here neutrino 
number profiles (r) with a finite virial radius rvir such 
that the total number of neutrinos is well defined 
and it is related to N^, and c according to Eq. For 
cosmological purposes we need to compute the average 
neutrino mass in the lump 

47r r^" 

("^1.) = ^ / drr'^n„{r)m,y{(l){r)). (26) 

The gravitational potential of the lump for large dis- 
tances r > Tvir is determined by (m^) and N^, as 



'M2 r 



(27) 



The overall size of the cosmologically averaged value of 
the fluctuations of the neutrino induced gravitational po- 
tential determines the size of possibly observable signals 
for neutrino lumps, as the integrated Sachs- Wolfe (ISW) 
effect in the CMB anisotropies [s^ or the bulk flow of 
peculiar velocities [i^ . 

In particular, we are interested in the average reduc- 
tion factor rm{Nv) = {fnv)/°nii,Q of the neutrino mass 
within lumps. We present in Fig. |S] the factor r„i as 
a function of the number of neutrinos in the lump. 
We show our results for two values of the concentra- 
tion parameter c = (10,30) to which corresponds iV* = 
(0.7, 0.4)iV^ and we assume a self similar relation between 
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lumps such that 

assuming 6 = 8.2 and b = 16.4. This relation arises for a 
fixed Us using Eqs. (fTTI) and (fT6|) . For a lump with iV* = 
4 X 10~^ and c = 30, these values of b lead, respectively, 
to the values of — l.SMpc and = 2.6Mpc that we 
used in previous sections. We find a substantial reduction 
{rjn < 1/2) even for relatively small lumps with ~ 
10~^A'tot. The spread of the different lines for the various 
sets of parameters can be interpreted as an indication 
of the uncertainty due to our limited knowledge of the 
actual profile n^{r). The qualitative behavior, however, 
is common to all cases. 

For lumps with number of neutrinos fraction Ni,/N tot 
in the range 10~^ to lO"""^ the average neutrino mass is be- 
tween one and three orders of magnitude smaller than its 
present cosmological value. This result demonstrates the 
significance of non-linear effects in the growing neutrino 
scenario, here quantified in terms of the backreaction on 
the neutrino mass. 

Our result implies important modifications for the 
gravitational potential of a neutrino lump when com- 
pared to linear perturbation theory. Since the neutrino 
mass is its source on all scales, the magnitude of the grav- 
itational potential is reduced by the same factor rm{Nu) 
compared to extrapolations based on the cosmological 
neutrino mass. 

By a similar argument, we expect an effective sup- 
pression of the cosmon mediated attractive force between 
neutrino lumps which is proportional to 2/3^. This can 
be estimated by the averaged interaction strength (/?) for 
the neutrinos in a neutrino lump, defined as 

The non-linear effects reduce the attractive force of the 
lump on a surrounding neutrino fiuid by a factor (/3) //3 as 
compared to the same number of neutrinos in the absence 
of non-linear effects. The attraction between two equal 
lumps is correspondingly reduced by a factor 
Furthermore, the characteristic time scale for the infall 
of neutrinos into a lump governed by the fifth force is 
inversely proportional to the interaction strength. This 
implies that the characteristic time scale for the infall is 
increased by a factor /3 / (/?) compared to the consider- 
ation excluding non-linear effects and thus results in a 
slow down of the infall. In the same line, the time scale 
for the clumping of lumps to larger ones is enhanced by 
a factor 

In Fig. [HI we show the averaged interaction strength 
(/3) as a function of the number of neutrinos in the lump 
for c — (10,30), again assuming b = 8.2 and b = 16.4. 
In the range N^/Ntot = 10~^ to 10~^, we find the aver- 
aged interaction strength (/3) to be reduced by up to one 
order of magnitude compared to the cosmological value 
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FIG. 8: Dependence of the ratio between the average neutrino 
mass in a lump and the background value on the number of 
neutrinos in a lump for a step NFW neutrino profile. 



/3 = —52 for a step NFW neutrino profile. Again, this 
suppression is a reflection of the non-linearities in the 
field equation. We found, however, that within the same 
radius R = erg, lumps with a pure NFW neutrino profile 
and with N^{R)/Ntot — 10~^ — 10~^ have an averaged 
interaction strength < 1, which implies that the ef- 
fective cosmon-neutrino interaction strength may even be 
smaller than the gravitational strength. 

In any case, if a large fraction of the neutrinos in the 
horizon is captured in large neutrino lumps, this may 
first alter the cosmological field equations for the back- 
ground fields by order one effects. Second, it may signif- 
icantly modify the growth of neutrino fluctuations with 
long wavelengths due to backreaction effects from fluc- 
tuations with shorter wavelengths. More precisely, ow- 
ing to neutrino clustering on smaller scales, fluctuations 
on large length scales will experience a modified effec- 
tive neutrino mass. In addition, they will be subjected 
to a smaller effective cosmon force and thus exhibit a 
substantially reduced growth rate compared to linear ex- 
trapolations. Our results for the gravitational potential 
can be used to compute the impact of neutrino lumps on 
the large scale structure and the cosmic neutrino back- 
ground in the form of the Integrated Sachs- Wolfe effect 

In summary, a fluid of neutrino lumps behaves very dif- 
ferently from a fluid of smoothly distributed neutrinos. 
If the typical size of the lump is sufficiently large and if 
the true neutrino number distribution in the lump should 
result in {(3) < 1, then the fiuid of lumps behaves quite 
similarly to dark matter, with only a moderate enhance- 
ment of the total attractive force. The time scale for 
a further growth of structure is comparable to the time 
scale for the gravitational growth of dark matter struc- 
tures. The early very rapid growth of neutrino structure 
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FIG. 9: The average interaction strength (/3) as a function 
of the number of neutrinos in the lump for a step NFW 
neutrino profile. 



terminates effectively once the lumps reach a critical size. 
A detailed estimate of this critical size would be highly 
interesting as it limits the size of the largest structures 
that can be expected in our universe. 

We also note that the relation between cosmological 
and terrestrial bounds for the neutrino mass may strongly 
depend on whether our galaxy is within a large neutrino 
lump or not. If our galaxy is indeed in a neutrino lump, 
then the choice of parameters (5 = —52 and rh in Eq. ^ 
(corresponding to a present cosmological neutrino mass 
fa^Q — 2.3eV) would lead to a substantially smaller neu- 
trino mass seen on a terrestrial experiment. Conversely, a 
present terrestrial bound < 2eV could be compatible 
with a substantially larger mass for unbound neutrinos, 
allowing realistic cosmologies with smaller values of 



V. THE VARIATION OF FUNDAMENTAL 
PARAMETERS 

A scalar field is expected to have couplings with other 
forms of matter (unless some unknown symmetry princi- 
ple explicitly forbids them) and consequently, if the field 
is varying, it generates variations of masses and coupling 
constants [11,11, 113, El]. In contrast to the neutrino- 
cosmon coupling, the interaction between the cosmon and 
baryons, leptons or photons must be much weaker than 
gravity in order to avoid confiict with the severe obser- 
vational bounds. Some suppression mechanism for the 
field dependence of particle masses (except for neutrinos 
and possibly dark matter) and dimensionless couplings 
as gauge or Yukawa couplings must be at work, perhaps 
by a fixed point behavior [ssl . l49j . Unfortunately, the 
strength of such a suppression is very model dependent 
such that a variation of a coupling of an observable size 



is not always guaranteed. However, if the neutrino mass 
arises from the field dependence of a heavy triplet in the 
cascade (or seesaw II) mechanism, this directly induces 
a small field dependence on the Fermi and fine structure 
constants ^22i] . Among the various fundamental couplings 
of the standard model of particle physics, we concentrate 
here on the fine structure constant a. Results for other 
couplings are similar. 

The possibility that dark energy couples to electro- 
magnetism and that it leads to the variation of the fine 
structure constant has been raised numerous times in 
the literature (13, E^. Claims that the fine structure 
constant was smaller in the past between redshifts 2-3, 

Aa/a 5 X 10"^ suggested by the Keck/HIRES high 

resolution quasar spectra, have been reported [50|, 
though these claims are under intense scrutiny and de- 
bate On the other hand, geological bounds from 
the Oklo natural reactor active at an equivalent redshift 
z « 0.14 yield Aa/a = (0.7 ± 1.8) x 10"^ [H and an 
independent analysis gives Aa/a — (0.6 ± 6.2) x 10~^ 
[5^ . For a linear dependence of Aa/a on redshift, we 
can immediately identify an inconsistency between the 
Oklo bounds and the claimed cosmological variation at 
high redshift. This discrepancy is approximately of a fac- 
tor 10^/(2/0.14) w 70. Additionally, strong constraints 
on the current variation arise from atomic clock measure- 
ments, a/a = (-1.6 ±2.3) x lO'^^yr"! Under the 
assumption that any observed variation of the fine struc- 
ture constant of the order 10~^ at redshifts 2-3 arises 
from a slowly time varying field which exhibits a small 
linear coupling to electromagnetism, the atomic clocks 
and the Oklo bounds challenge, therefore, models of dark 
energy with a monotonic scalar field evolution f56i-60] . 

There are, however, several factors that may allevi- 
ate the apparent discrepancy between the bounds on the 
present variation and the Oklo bounds and the claims of 
a variation at high redshift. For example, the time evolu- 
tion of the field may have slowed down in recent times due 
to the coupling of the quintessence field to matter [6l|, 
such as in thegrowing neutrino scenario or in crossover 
quintessence [62l- We discuss here another effect, namely, 
that our galaxy may be within a large neutrino lump 
such that the cosmon field has decoupled from the cos- 
mological evolution at a recent redshift, say z w 1. We 
expect then a substantial reduction of the present time 
variation and predict a different value of the fundamen- 
tal constants in different regions of the lump and within 
different lumps. Indeed, we will see that inside a lump 
the time variation of the cosmon field can be sufficiently 
small to prevent a violation of the Oklo bound even for 
the substantial variation of a claimed in Refs. [sol [5l|. 

In the following, we assume that we live within a large 
neutrino lump which formed rapidly around z ~ 1 — 2 
and since then it maintained an approximately constant 
neutrino profile as given in Eq. ([9]) or ((T5|) . We will iden- 
tify a set of model parameters for the growing neutrino 
scenario, which effectively leads to a decoupling of the 
value of the scalar field inside a neutrino lump and its 
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cosmological value. 

In the case in which the evolution of the scalar field is 
sufficiently small, i.e. |A0| = |0 — 0o| < M, we can take 
the cosmon coupling with electromagnetism to be linear. 
This can be seen as the first term of a Taylor expansion. 
The evolution of a is then given by 



Aa _ a — ao 
a ao M 



(30) 



where C is a constant. 

Accordingly, the local variation of a (by local we mean 
as measured with atomic clocks or with Oklo) is given by 
the difference of the local value of the field at two different 
epochs (2: w or z = 0.14 for atomic clocks and Oklo, 
respectively). Hence, for Oklo, we have 



Aa 

a 



Oklo 



M 



0.14, r) - 0(z = 0,r)] , (31) 



where r represents our radial position in the neutrino 
lump. Similarly, the cosmological variation of a can be 
seen as a result of the difference between the value of the 
field at a distance Rc from the center of the lump (where 
the observed molecular cloud is located) at high redshift 
(e.g. Zc = 2) before the lump has formed and the local 
value of the field today 



Aa 

a 



cloud 



M 



[cf>{zc = 2,R,)-(l)iz = 0,r)] . (32) 



For practical purposes we will first identify 4'{zc, Rc) with 
the cosmological value ^(zc) a-nd we will discuss other 
scenarios for (l){zc, Rc) below. 

In order to solve the discrepancy between the Oklo 
bound and the claimed variation of a cosmologically, the 
ratio between Eqs. (|3T|) and ([32]) 



/(r) = '^(! = 0-14,r)-(/)(z = 0,r) ^ ^^^^ 



(ze = 2)-0(z = O,r) ' 



must satisfy 



4.4 X lO"-' < /(r) < 1.9 X 10"% 



(34) 



where we have used the strongest Oklo bound of Ref. [53| 
and assume that Aa/a{zc = 2) = —5.7 x 10"^ as sug- 
gested by the results of Refs. [sollsH. 

The quantity / above can be thought of as the ratio of 
two distinct effects. On the one hand we have a reduction 
factor between the local evolution of the field (between 
redshift zi of the local observation and today) and the 
background evolution 



^ 4>{zi,r) - 4){z = 0,r) 
4>{zi)-4>{z = Q) 



(35) 



This quantity is, up to the normalization by the back- 
ground, identical to A</) illustrated in Fig. 01 On the 
other hand we have a measurement of the difference be- 
tween the value of the field outside the lump in the distant 
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FIG. 10: The suppression factor /(r) for four neutrino num- 
ber profiles with TV, = 4 x 10~^. The horizontal dashed lines 
represents the values between which the claimed observed cos- 
mological variation of a and the Oklo bounds are consistent. 

past at redshift Zc and its present value compared to the 
cosmological evolution in the same time interval 

/.(^c,.,i?.)^%^f^m^. (36) 

(t>[Zc) - (p(z = 0) 

Then we obtain that for a generic zi and Zc, the ratio / 
is simply given by 



fi{zi,r) ${zi)-^{z = 0) 



fc{zc,r,Rc) (j){zc) - (j){z ^ Q) 



(37) 



Provided the value of the suppression factor // is suffi- 
ciently smaller than the cosmological shift factor /c, local 
bounds at zi become easier to be satisfied even for a large 
cosmological variation of a at Zc. We emphasize that the 
same ratio /(r) applies to all time variations of couplings 
to which a linear approximation similar to Eq. pO|) holds. 

In Fig.[Tn]we show the overall suppression factor /(r) 
in Eq. ([55)1 for our model. We also indicate with the 
dashed lines the value of / between which the discrep- 
ancy between the Oklo bounds and claimed cosmological 
variation is lifted. In our example with A^* = 4 x 10""^, 
no discrepancy is present if our local position is r < lOrs 
for the pure NFW profile whereas for the step NFW, r 
is constrained to be r < 3rs. 

The suppression factor depends strongly on the num- 
ber of neutrinos in the lump. In Fig. [TT] we show 
/(r = O.Its) as a function of N^/Ntot for b — 8.2 and 
c = 30. It becomes evident that only for large lumps 
the claimed observed cosmological variation of a and the 
Oklo bounds is consistent. 

Finally, we notice that the field evolution outside the 
lump, (/'(z) — 0(0), is a useful quantity for comparison, 
but actually does not appear in the observable quanti- 
ties. For the local observations, including Oklo, only the 
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FIG. 11: The suppression factor /(r = O.Irs) for the step 
NFW profile with respect to the number of neutrinos in a 
lump for b = 8.2 and c = 30. The horizontal dashed lines 
represent the current Oklo bounds. 



difference (j){zi, r) — 0(0, r) matters, while the observable 
cosmological variations only depend on (^(zc i?c)— 0(0, r). 
In this context we note that (f>{zc, Re) depends on the ob- 
served position of the molecular cloud. If lumps form 
early enough, it is even conceivable that at a cosmo- 
logical redshift Zc ~ 1, the scalar field is already suffi- 
ciently inhomogeneous such that a possible cosmological 
time evolution of couplings is replaced by a dependence 
of Aa/a on the position of the molecular cloud. This 
would result in a dependence of Aa/a on the particular 
observed quasar. In the case of an inhomogeneous sit- 
uation, 0(zc,-Rc) ought to be replaced by (f>{zc, Rc,0,v), 
i.e., it depends on the angular coordinates. (For large 
enough Zc the distance Rc to the cloud can be related to 
the redshift of the cloud Zc-) A systematic survey of the 
variation of fundamental couplings with respect to these 
coordinates would produce a map of the inhomogeneities 
of the cosmon field. An analysis suggesting an angular 
dependence of a has recently been reported [631] . 



VI. CONCLUSIONS 

In this work, we have explored backreaction effects 
arising in growing neutrino quintessence. They result 
from the response of local scalar field variations to the 
formation of large scale non-linear neutrino lumps in the 
presence of a cosmon- neutrino coupling. The coupling 
strength is typically much larger than the gravitational 
interaction. This amplifies the effects of the non-linear 
nature of the neutrino source term in the field equation 
inside of a highly-concentrated neutrino lump. By nu- 
merically solving for the spatial scalar field profile as- 
suming an NFW distribution for the neutrino number 



density, we have demonstrated that these non-linearities 
are manifested in a strong reduction of the neutrino mass 
within a neutrino lump. Stunningly, for realistic exam- 
ple cases, we found in the interior of a lump, a suppres- 
sion of the neutrino mass by several orders of magnitude 
compared to the cosmological value well outside of the 
lump. A key result of our work is that the radial profile 
of the cosmon field and thus of the neutrino mass are 
effectively frozen inside a lump, once it virializes, while 
both the cosmological field and the background neutrino 
mass may continue to evolve with time. The radial de- 
pendence of the scalar field and the neutrino mass in the 
interior of the lump are essentially fixed by the neutrino 
distribution. 

For the time being, substantial uncertainties regarding 
the neutrino number distribution in the lump subsist. 
We have investigated here an NFW profile, but it should 
be mentioned that neutrino lumps differ from dark mat- 
ter clusters in several important aspects: (i) the relevant 
time scale of formation is much shorter, such that there 
is plenty of time to reach a quasi-static equilibrium con- 
figuration; (ii) the varying mass of the neutrinos can have 
important effects for the stable profile; (iii) if the neutrino 
mass in the core of the lump is very small, the neutrinos 
with high enough momentum behave close to relativistic 
particles, leading to a less concentrated profile. The true 
neutrino profile may be determined by a self-consistent 
stable static solution similar to Refs. [63 - l66j . For this 
purpose, however, one needs to gain some insight on the 
average velocity distribution of the neutrinos inside the 
lump, resulting in an effective pressure which depends on 
the radius. The strong reduction of the neutrino mass in- 
side the lump is an effect that seems to be rather robust, 
independent of the details of the neutrino number dis- 
tribution. The size of the effect mainly depends on the 
total number of neutrinos contained in the lump, and 
to a smaller extent on the core density or concentration. 
We have demonstrated this by investigating distributions 
with different neutrino numbers and core densities. In 
principle, also the neutrino number distribution can be 
influenced by a change of the cosmological evolution of ip. 
We have neglected this effect taking the neutrino num- 
ber distribution to be static after virialization. In view 
of the effective decoupling of the lump from the cosmo- 
logical evolution this ansatz is self-consistent, at least as 
far as the qualitative features are concerned. 

We have computed the average neutrino mass {m^) 
within various families of neutrino lumps. For the specific 
scenario under consideration and for lumps with N^/Ntot 
in the range 10^^ to 10~^, we found its value to be be- 
tween one and three orders of magnitude smaller than the 
cosmological mass. This result reflects the importance of 
backreaction effects in growing neutrino quintessence. In 
particular, it causes the gravitational potential of very 
large neutrino lumps to also be reduced by several or- 
ders of magnitude when compared to estimates based on 
the value of the cosmological neutrino mass. This reduc- 
tion stems from the fact that the gravitational potential 
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is sourced by the total neutrino mass in the lump and 
therefore sensitive to the average neutrino mass. Sim- 
ilarly, we found a substantial reduction of the effective 
coupling between large neutrino lumps and the cosmon. 

Finally, we have considered the implications of our 
findings for the variation of the fine structure constant 
a, which would result from a coupling of the cosmon to 
electromagnetism, or the variation of other fundamen- 
tal constants. Variations of particle masses and coupling 
constants result from variations of the value of the cos- 
mon field (f>. We find that the time variation of (f) is sub- 
stantially reduced inside large neutrino lumps, as com- 
pared to the outside variation. If our galaxy is located 
within such a lump, local time variations - as observed by 
precision experiments or the analysis of the Oklo natural 
reactor - would be substantially suppressed as compared 
to cosmological observations. This may weaken the dis- 
crepancy between bounds for local variations of a and 
claims of an observed variation in molecular absorption 
spectra at high redshifts, which would result from sim- 
ple linear extrapolations. For sufficiently strong inhomo- 
geneities of the cosmon field at a redshift around z — 2, 
the observed values of Aa/a could even depend on the 
position of the molecular clouds and therefore on direc- 
tion. In our model, the maximal variation of a is limited 
by the maximal ratio of the neutrino masses outside and 
inside our lump. The neutrino mass locally must exceed 
nil, > 0.05eV and its value in the background cannot 
be larger than a few eV without confiicting with cosmo- 
logical limits. These bounds consequently constrain our 
position inside the lump. 

Our investigation illustrates the significant impact of 
strong backreaction effects in the growing neutrino sce- 
nario on the interpretation of cosmological observations. 
These non-linear effects are expected to be even larger in 
the case of a field dependent cosmon-neutrino coupling 
/3(0). Unfortunately, the large size of the backreaction 
effects also implies that only once backreaction is incor- 
porated properly in the full analysis, the range of param- 
eters of growing neutrino quintessence, which is compat- 
ible with cosmological observations, can be analyzed. In 
particular, the growth of large scale neutrino perturba- 
tions itself can be expected to be significantly modified. 
The reason is that after the formation of neutrino lumps 
on smaller scales these large scale perturbations feel a 
reduced effective neutrino mass and effective coupling. 

Our findings have important consequences for the ef- 
fects of neutrino lumps on the cosmic microwave back- 
ground in the form of the Integrated Sachs Wolf effect 
[38j or on enhanced bulk flows [46j. Indeed, neutrino 
lumps would be observable mainly by the gravitational 
potential induced around them. In this sense, they are 
very similar to dark matter clusters. Galaxies will have 
tendency to fall into the large scale lump potentials. The 
time when this effect becomes important and the size of 
the effect depends crucially on the size of the gravita- 
tional potential generated by the lumps and therefore on 
the effective neutrino mass inside the lump. Also the cos- 



mological background evolution will be modified once a 
large fraction of the neutrinos is within lumps. 

The growing neutrino scenario seems to be one of 
the first examples for a realistic backreaction effect 
which modifies the parameters of the homogeneous and 
isotropic field equation by order one effects. 



Appendix A 



In this appendix, we consider neutrino number density 
profiles Hi, (r) with a finite virial radius rvir such that the 
total number of neutrinos iV^ is well defined. In this case 
we can relate the parameters used in the main text to 
the neutrino number overdensity A„ = {n^)/n^o. 

The neutrino structure, bound by the fifth force and 
by gravity [23l l63 - [66| has an average density (n^) defined 
by 



An 



(Al) 



We can rewrite this equation in terms of the ratio of the 
number of neutrinos in the lump, iV^, with respect to the 
total number of neutrinos, A^tot, within the horizon. 



1 



Hq V Atot A„ 



1/3 



(A2) 



where A'tot — 47rniyo/(3iJo)- The conservation of the 
number of neutrinos taking part in the formation process 
of the lump and the self-similarity of the free-fall of col- 
lapsing structures may suggest that A„ is a constant. In 
this case, however, the ratio of the average energy density 
{pu) of the lump to the present neutrino energy density in 
general does not coincide with A„, {pv)lpvO ^ A„. The 
reason is that with increasing density of the forming lump 
the mass of the neutrinos in the interior increases slower 
than the background neutrino mass or it can even freeze 
or decrease. If A„ is a universal constant for all types 
and sizes of lumps, then it partly fixes the parameters 
Ts and Us of the neutrino number density distribution. 
Inserting Eq. (|A2[) into Eq. we can determine the 
functional form of the corresponding NEW distribution 
for a given Ni, in terms of A„ and the concentration c. 



(ci?o) 



-M.otA, 



47ri^(c)' 
1 / A, 



cHo V A^totA„ 



1/3 



(A3) 
(A4) 



In Eig. [T2]we plot the resulting functional dependence 
of the virial radius rvir on the number of neutrinos N^, 
it contains for A„ = 1.1 x 10"* according to Eq. (jA2p . 
For the dedicated numerical simulation in Sec. IIIIl the 
neutrinos contained in a volume with radius R = 25Mpc 
have reached this overdensity at z — 2.017. The spherical 
infall of the neutrino shells corresponding to radii r < R 
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FIG. 12: The virial radius of a neutrino lump as a function 
of the total number A'^^ of neutrinos it contains for A„ — 
1.1 X 10". 



at z = 2.017 is not affected by neutrinos in shells with 
r > R. We may thus discuss a situation with no neutrinos 
for r > R at z — 2.017. The dynamics of the neutrinos 
in the shells with r < i? is the same as in the simulation. 
For this scenario the value A„ = 1.1 x 10^ constitutes 
a lower bound for the final overdensity as A„ can only 
increase due to a still ongoing collapsing process. 
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